Draft version January 4, 2012 

Preprint typeset using L^-T^X style cmulateapj v. 5/2/11 



THE EFFECT OF DIFFERENT OBSERVATIONAL DATA ON THE CONSTRAINTS OF COSMOLOGICAL 

PARAMETERS 

Yungui Gong 1 ' 2,3 , Qing Gao 2 , and Zong-Hong Zhu 4 

Draft version January 4, 2012 

ABSTRACT 

The constraints on ACDM model from type la supernova data alone and BAO data alone are 
similar, so it is worthwhile to compare their constraints on the property of dark energy. We use the 
SNLS3 compilation of 472 type la supernova data, the Gamma Ray Bursts data, the baryon acoustic 
oscillation measurement of distance, the cosmic microwave background radiation data from the seven 
year Wilkinson Microwave Anisotropy Probe, and the Hubble parameter data to study the effect of 
their different combinations on the fittings of cosmological parameters. Neither BAO nor WMAP7 
data alone gives good constraint on the equation of state parameter of dark energy, but both WMAP7 
data and BAO data help type la supernova data break the degeneracies among the model parameters, 
hence tighten the constraint on the variation of equation of state parameter w a , and WMAP7 data 
does the job a little better. Although BAO and WMAP7 data provide reasonably good constraints 
on O m and O^, it is not able to constrain the dynamics of dark energy, we need SNe la data to probe 
the property of dark energy, especially the variation of the equation of state parameter of dark energy. 
The addition of H(z) data helps better constrain the geometry of the universe fifc and the property 
of dark energy. For the SNLS SNe la data, the nuisance parameters a and /3 are consistent for all 
different combinations of the above data. Their impacts on the fittings of cosmological parameters 
are minimal. ACDM model is consistent with all the observational data and it is favored against 
Dvali-Gabadadze-Porrati model. 

Subject headings: cosmological parameters; dark energy 



1. INTRODUCTION 

The accelerating expansion of the universe was first 
discovered in 1998 by the observatio n s of Type la su- 
pcrnovac (SNe la) (|Riess et al.l 119981: iPerlmutter et al.l 
1999). As more accurate data are available, it is possible 
to measure the acceleration and the dynamical mecha- 
nism behind the acceleration. There are three different 
possibilities for the acceleration. The first possibility is 
that a new exotic form of matter with negative pressure, 
dubbed as dark energy drives the Universe to acceler- 
ate. The cosmological constant is the simplest candi- 
date of dark energy which is also consistent with obser- 
vations, but at odds with quantum field theory. The 
second possibility is that general relativity is modified at 
the cosmologic al scale, such as Dvali-Gabadadze-P orrati 
(DGP) model (jDvali. Gabadadze fc PorTati|[2000t) . The 
third possibility is that the universe is inhomogeneous. 
In this paper, we consider the possibility of dark energy 
only. 

In the recent release of the measurements of the baryon 
acoustic oscillation (BAO) peaks at redshifts z — 0.44, 
0.6 and 0.73 in the galaxy correlation functio n of the final 
dataset of the WiggleZ dark energy survey, iBlake et"aTl 
(2011) used these three BAO data along with BAO 
data at redshifts z — 0.2 and 0.35 measured from the 
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distribution of galaxies (jPercival et al.l 120101 ) and the 
measurement of BAO at redshift z = 0. 106 from the 
6-deg ree Field Galaxy Survey (6dFGS) (jBeutler et al.l 
1201 If) to constrain ACDM model. It was found that 
the constraints on ACDM model from BAO data only 
are even better than t hose from Union2 SNe la data 
(jAmanullah et al.l 120101 ) only. This means that the cur- 
rent BAO data is robust to constrain cosmological pa- 
rameters. Because there are more than 500 SNe la data 
points, and only 6 data points in the updated BAO data, 
the fitting process for BAO data is easier and faster than 
that for SNe la data. If the constraints on the property 
of dark energy from BAO data are tighter than those 
from SNe la data, then we just need to apply BAO data 
only for a faster fitting although SNe la data and BAO 
data are complementary to each other. In this paper, we 
use a simple dark energy model to test the robustness of 
BAO data, and compare the constraints on the equation 
of state of dark energy from different data. 

The geometry of the universe is sensitive to the cos- 
mic microwave background (CMB) data, so it is ex- 
pected that the addition of the seven-year Wilkinson 
Microwave Anisotropy Probe (WMAP7) data to cither 
SNe la or BA O data helps tigh ten the constraints on 
ACDM model. IBlake et al.1 (|2011[) found that the combi- 
nation of BAO and WMAP7 data gives much better con- 
straints on ACDM model t han the combination of SNe 
la and WMAP7 data does. ISullivan et all (|2011D found 
that both SNLS3 SNe la data alone and the combina- 
tion of WMAP7 and BAO da ta at redshifts z = 0.2 and 
z = 0.35 (|Percival et al.lr2~010D gave similar constraint on 
the equation of state parameter w for the flat constant w 
model. The redshifts of BAO data span from z — 0.106 
to z = 0.73, we may expect that BAO data catches the 
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dynamical property of dark energy, so it is necessary to 
study the effects of different observational data and their 
combinations on the constraints on the equation of state 
of dark energy. 

The question whether dark energy is just the cos- 
mological constant remains to be answered. Re- 
cently, there are lots of studies in determining 
wheth er ACDM model is consistent with observa- 
tions (iHuang et al.|[ 2009; Shaficl oo. Sahni fe Starobinskyl 
2001 ICai Su fc Zhand 120101: iLampeitl et al.l 1200a 
Serra et alJl2009b jGong et al.ll2010aHGong. Wang fc Cail 
2010bHPan eTaLll2010l : IGong Zhu fc ZhiJl20lH ILieTaLl 
201 If) . iLi. Wu fc"Yul d2011|) considered the tensions 
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between different dataset through the reconstruction 
of Om(z) by usi ng Chevallier-Polarsk i- Linde r (CPL) 
param etrization (Chcvallie r fc Polarskil 120011 : iLinderl 
2003) of the equation of state of dark energy. They found 
that a tensio n between low r edshift and high redshift 
data existed. ICai. Su fc Tuol (|201lD used the fi gure of 
merit (FOM) propos ed by the Dark Energy Task Force 
(|Albrecht et al.|[200l) as a diagnostic to study the effec- 
tiveness of different combinations of data on constraining 
wq and w a in CPL model. 

In this paper, we first study the robustness of BAO 
data, then study the constraints on the equation of state 
of dark energy based on different combinations of the 
following data: the three year Supernova Legacy Sur- 
vey (SNLS3 ) sample of 472 SN e la data with system- 
atic errors (jConlev et al.l 120111): the 59 Gamma Ray 
Bursts fGRB) data (iWeil 120101): the BAO measure- 
ments from the 6dF GS dBeutler et al.l 1201 ID . the dis- 
tribution of galaxies dPercival et al.l | 2010l) and the Wig- 
gleZ dark energy survey (iBlake et al.boTll ): the WMAP7 
data (Ko matsu et al.l 1201 If) : and the Hubble parameter 
H(z) data (jGaztahaga. Cabre fc Huill2b09bl : iStern et al.l 
2010). In addition to studying the effects of different 
observational data and their combinations on the con- 
straints of cosmological parameters, we also reconstruct 
the equation of state of dark energy w(z), the decelera- 
tion parameter q(z) and Om{z) by using these datasets. 

The paper is organized as follows. In section 2, 
we p resent the SNLS3 SNe la d ata (jConlev et al.1 
f201lh. the GRB d a ta dWell l2010l), the BAO data 
dBeutler et ail 120111: iBlake et al.l 120111: IPercival et al.l 
|2010|), the WMAP7 data (iKomatsu et al.l 1201 ID. the 
H(z) data (jGaztahaga. Cabre fc Huill2009bl IStern et al.l 
I2010T) . and all the formulae related to these data. In sec- 
tion 3, we present all the models and the fitting results, 
and conclusions are drawn in section 4. 



2. OBSERVATIONAL DATA 

The SNLS3 SNe la data consists of 123 low-redshift 
SNe la data with z < 0.1 mainly from Calan/Tololo, 
CfAI, CfAII, CfAIII and CSP, 242 SNe la over the red- 
shift range 0.08 < z < 1.06 observed from the SNLS 
(jConlev et al.l 1201 ID . 93 intermediate- redshift SNe la 
data with 0.06 < z < 0.4 observed during the first sea- 
son of Sloan Digital Sky Survey (SDSS)-II supernova 
(SN) survey (jKessler et al.l l2010h . and 14 high-redshift 
SNe la data with z > 0.8 from Hubble Space Tele- 
scope (jRiess et al.ll2007h . The SNLS3 SNe la data used 
the combination of SALT 2 and SiFTO light-curve fitters 
(jConlev et alJ[20H . To use the 472 SNLS3 SNe la data 
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C 1 ) 

where tub is the rest-frame peak B-band magnitude of a 
SN, the predicted magnitude of the SN given a cosmo- 
logical model is m mod = 5 log 10 V L (z heU z cmb , p) - a(s - 
1) +/3C + Mb, Zhei and z cm b are the heliocentric and the 
CMB frame redshifts of the SN, s is the stretch given 
by the data, C is the color measure for the SN given by 
the data, a and /3 are nuisance parameters used for the 
SNLS3 data fitting, Mb is another nuisance parameter 
incorporating the absolute magnitude and Hubble con- 
stant and it is marginalized over in the SN fitting process 
because of the arbitrary normalization of the magnitude, 
C sn (zi, Zj) is the covariant matrix which includes both 
the syste matical and statist ical uncertainties for the SNe 
la data (jConlev et al.l [20TlT ) . The correction on the de- 
pendence of the host-galaxy stellar mass is also included. 
The Hubble-constant free luminosity distance T>l(z) is 



V L (z) = H d L (z) = ±±f= S k 
V 



dx 



E(x) 



> (2) 



where the dimensionless Hubble parameter E{z) = 
H(z)/Ho, and Sk(x) is defined as x, sin(x) or sinh(x) 
for fc = 0, +1, or -1, respectively. For the fitting to the 
SNL3 data, we need to add two more nuisance parame- 
ters a and f3 in addition to the model parameters p and 

the nuisance parameter Mb- 

The fitting of GRB data (jWeill2010f ) is similar to that 
of SNe la data except that the distance modulus ^ = 
m — M = 51og 10 T>l(z, p) is used instead, so the fitting- 
is simpler. 

For t he BAO data, we us e the measurements from the 
6dFGS dBeutler et~al1l201lD . the distribution of galaxies 
(IPercival et al.ll2010D a n d the WiggleZ dark energy sur- 
vev (jBlake et al.ll201lD . IPercival et all (j2010l ) measured 
the distance ratio, 



r s {zd) 
D v (z) 
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at two redshifts z = 0.2 and z — u.jo w uc u, 2 
0.1905 ± 0.0061, and dg 6 | 5 = 0.1097 ± 0.0036. Here the 
effective distance is 
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z d is the drag redshift defined in Eis enstein fc Hul (jl998D . 
the comoving sound horizon is 



r s (z) 



c s (x)dx 



(5) 



where the sound speed c s (z) = l/-\/3[l + i?f,/(l + z)], 
and R b = m b h 2 /{A x 2.469 x 10" 5 ). iBeutler et all ([Ml 
derived that d%% 6 = 0.336 ± 0.015. The WiggleZ dark 
energy survey measured the acoustic parameter 
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at three redshifts z — 0.44, z — 0.6 and z = 0.73, and the 
results and t heir covarian c e mat rix are listed in table 3 
and table 2 in lBlake etabl (|2011ft . To use the BAO data, 
we minimize 

2 

xL ( P , ^th 2 , h) = MiC^idk, d 3 )Ad, 

(do.ioe - 0.336) 2 
0.015 2 

3 

+ J2 AAiC^iAuA^AAj, (7) 

where ck = (d z=0 .2, ^=0.35), = — d° 6s and the 
covariance matrix Cd?.(dj,d j) for at = (0.2, 0.35) is 
taken from equation (5) in IPercival et al.l (|2010l ); At = 
(A(0.44),A(0.6),A(0.73)), AA4 = A{ Zi ) - A{ Zi ) obs and 
the covariance matrix CA(Ai, Aj) for the data points 
A(z) at z = (0.4 4,0.6,0.73) is taken from table 2 in 
iBlake et"aT] ()2011|) . Besides the model parameters p, we 
need to add two more nuisance parameters fifth 2 and 
fl m h 2 when we use the BAO data. 

For the WMAP7 data, we use the measurements of the 
three derived quantities: the shift parameter R(z*) and 
the acoustic index Ia(z*) at the recombination redshift 
z* . In particular, we minimize 

3 

XcMB(P^bh 2 ,h) = AxiCcM B (xi,Xj)Axj, (8) 

where the three parameters Xi = [R(z*), Ia(z*), z*], 
Axi = Xi — x° bs and the covariance matrix Ccm b (%i , %j ) 
for the three p arame ters is taken from Table 10 in 
IKomatsu et al.l (|2011D . The shift parameter R is ex- 
pressed as 



R(z*) 



^UV L {z*) 



1 + z* 

The acoustic index I a is 
ird L (z*) 



l A {z*) 



1.725 ±0.018. 



302.09 ±0.76, 



(9) 



(10) 



(1 + z*)r s (z*) 

and z* is the redshift at the recombinati o n wit h the 
parametrization defined in IHu k Sugivamal (|1996f ). We 
also need to add the nuisance parameters fifth 2 and fl m h 2 
to the parameter space when we employ the WMAP7 
data. 

Additionally, we use the H(z) data at 11 different red- 
shifts obtai ned from the differ ential ages of red-envelope 
galaxies in iStern et "all (|2010f ) , and three more Hubble 
parameter data H(z = 0.24) = 76.69 ± 2.32, H(z = 
0.34) = 83.8 ± 2.96 and H(z = 0.43) = 86.45 ± 3.27, 
determined by IGaztanaga. Cabre k H ui (20093). So we 
add these H(z) data to % 2 , 



X 2 H (P,h) 
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H bs(zi)Y 
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where a hi is the la uncertainty of H(z). Basically, the 
model parameters p are determined by minimizing 

(12) 



The likelihood for the parameters p in the model and 
the nuisance parameters is computed using the Monte 
Carlo Markov Chain (MCMC) method. The MCMC 
method randomly chooses values for the above param- 
eters p, evaluates \ 2 an( i determines whether to accept 
or reject the set of parameters p using the Metropolis- 
Hastings algorithm. The set of parameters that are ac- 
cepted to the chain forms a new starting point for the 
next process, and the process is repeated for a suffi- 
cient number of steps until the required convergence 
is reached. Our MCMC code i s based on the pub- 
licly available package co smomc (|Lewis & Bridle! 2002: 
IGong. Wu k Wandl2008l) . 

After fitting the observational data to different 
dark energy models, we apply th e Om diagnostic 
(|Sahni. Shafieloo k Starobinskvil2008( ) to detect the de- 
viation from the ACDM model. F or a flat universe 
(jSahni. Shafieloo k Starobinsky|[200l . 
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For the flat ACDM model, Om(z) = Q m is a constant 
which is independent of the value of VL m . Because of 
this property, Om diagnostic is less sensitive to observa- 
tional errors than the equation of state parameter w(z) 
does. On the other hand, the bigger the value of Om(z), 
the bigger the value of w(z), so the behavior of Om(z) 
catches the dynamical property of w(z). 

We also apply the FOM as a diagnostic tool to com- 
pare the effectiveness of different combinations of obser- 
vational data on constraining the equation of state pa- 
rameters wq and w a in CPL model. FOM is defined as 
the the reciprocal of the area of the error ellipse enclos- 
ing the 95% confidence limit in the wo~w a plane, it is 
proportional to [detC w (wo,w a )}^ 1 ^ 2 , here C w (wo,w a ) is 
the correlation matrix of wq and w a . 

3. COSMOLOGICAL FITTING RESULTS 

We first review the effects of different combinations of 
observational data on ACDM model. The f2 m -fiA con- 
tour for applyi ng only the SNLS3 SNe data was shown 
in Figure 8 in iConlev et al.l (|2011l). By combining the 
SNLS3 SNe and the WMAP7 data. ISullivan et al.l (l20ll 
obtained the constraint on VL m and fife, and fi m -f2fc con- 
tour was shown in Figure 4 of their paper. From these 
results, we see that SNLS3 SNe data alone does not pro- 
vide tight constraint on J7 m and f^-. With the addition 
of WMAP7 data, the constraint on J\. becomes much 
tighter, hencefore tightens the constraint on f2 m . There- 
fore, WMAP7 data can be used to tighten the constraint 
on the geomet r y of t he universe as shown i n Figure 15 
in IBlake et al.1 (|2011| ). In the same figure, IBlake et al.1 
(2011) showed the constraint on f2 m and ilk by ap- 
plying the BAO data only with the assumption that 
fli>h 2 = 0.02227. Comparing the constraints from SNLS3 
or Union2 SNe la data alone with that from BAO data 
alone, we see that the constraints are similar, and the 
constraint on Vl m is much better by using BAO data 
alone than that by using SNe la data alone. IBlake et al.l 
(|2011f) also compared the constraints on ACDM model 
by using the combination of BAO and Union2 SNe la 
data with the addition of WMAP7 data, and they found 
that fl m -flk contour became much smaller with the com- 
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bination of BAO and WMAP7 data compared with that 
using the combination of Union2 SNe la and WMAP7 
data. Moreover, the constraint from the combination of 
Union2 SNe la, BAO and WMAP7 data is similar to 
that from BAO and WMAP7 data. These results show 
that BAO data mainly tightens the constraint on VL m 
and WMAP7 data mainly tightens the constraint on Q^, 
while current SNe la data still gave large f2 m -f2fc contour. 
For comparison, we show all the constraints in Figure 
[TJ Note that we set fi^/i 2 as a free nuisance parameter. 
In Figure HJ we show the constraints on ACDM model 
from SNLS3 SNe la data alone (the green lines), BAO 
data alone (the yellow line), the combination of SNLS3 
SNe la and BAO data (the cyan lines), the combina- 
tion of SNLS3 SNe la and WMAP7 data (the magenta 
lines), the combination of BAO and WMAP7 data (the 
blue lines) , the combination of SNLS3 SNe la, BAO and 
WMAP7 data (the red lines), the combination of SNLS3 
SNe la, BAO, WMAP7 and H(z) data (the black lines), 
and the combination of SNLS3 SNe la, BAO, WMAP7, 
GRB and H(z) data (the shaded regions). As expected, 
GRB data has little effect on the Q m -flk contour when 
it is co mbined with the SNe la data. This point was also 
found in lCai. Su fc Tuol (|2011| ). In the following analysis, 
we do not consider the effect of GRB data alone on con- 
straining cosmological parameters. H(z) data further re- 
duces the errors on Qk and moves the best fit value of f2fe 
toward zero. These results along with the constraints on 
the nuisance parameters a and /3 are summarized in ta- 
ble 1. The error bars of a and j3 are around 0.1, and they 
are consistent at la level for different fittings. The best 
fit value of £l m from SNe la data is marginally consistent 
with that from BAO data, this shows a little tension be- 
tween SN la and BAO data. The tension between BAO 
measurement and higher redshift type la supernova (SN 
la) was noticed inlPercival et all (l2007f) . and the tension 
was lessened in lPercival et al.l ((20100 due to revised error 
analysis, different methodology adopted and more data. 

3.1. CPL parametrization 

In this section, we apply the CPL p arametrization 
(|Chevallier fc Polarskil[200lt ILinderl 120031) . 

7/1 7 

w(z) =w + (14) 
1 + z 

to test the effects of different combinations of data on 
constraining the property of dark energy. In this model, 
we have four model parameters p = (fi m , £lk, wq, w a ). 
From the results of ACDM model, we know the con- 
straints on fifc from either SNe la data alone or BAO 
data alone are not good. For the curved CPL model, 
due to the addition of two more model parameters, we 
expect the situation becomes worse and in fact this is 
true. So we only consider the combinations of SNLS3 
SNe la and/or BAO data with WMAP7 data. The con- 
tours are shown in figure [2] and the upper panels of figure 
03 The lcr results along with the constraints on the nui- 
sance parameters a and (3 are summarized in table 2. 
As expected, we see that the constraint on flk from the 
combination of BAO and WAMP7 data (the blues lines) 
is better than that from the combination of SNe la and 
WMAP7 data (the magenta lines). The f2 m -fifc contour 
becomes much smaller when we combine SNe la, BAO 
and WMAP7 data (the red lines). The addition of H(z) 



further reduces the errors on flk and moves the best fit 
value of flk toward zero. From the wo~w a contours in Fig- 
ure [2] we see that the constraints from the combination 
of SNe la and WMAP7 data are much better than those 
from the combination of BAO and WMAP7 data, and the 
FOM is almost 5 times larger. The uncertainties in w a 
from the combination of SNe la, BAO and WMAP7 data 
are reduced more than half compared with those from the 
combination of SNe la and WMAP7 data, the FOM in- 
crease more than 4 times. When SNe la data is combined 
with other data, fife is degenerated with wq and w a . As 
flk increases, the uncertainties in w a reduced. However, 
for the combination of BAO and WMAP7 data, the de- 
generacies between flk and wq and w a become weaker. 
Therefore, although SNe la data does not provide tight 
constraints on Q m and Ofc, its constraint on the equa- 
tion of state parameter of dark energy is much better. 
ACDM model (the cross) is outside the la contour when 
we use the combination of SNe la, BAO and WMAP7 
data (the red lines). With the addition of H(z) data, 
the WQ-Wa contour is further reduced and ACDM model 
is inside the la contour. By using the constraints from 
the combination of all observational data, we reconstruct 
w(z) and the result is shown in figure 2] 

Since WMAP7 data greatly reduces the error on fifc, 
its effect on flat CPL model still needs to be studied. 
In Figure [H we show the marginalized wo-w a contour 
plots constrained from different combinations of data. 
The la results along with the constraints on the nui- 
sance parameters a and {$ are summarized in table 2. 
The constraints on wq and w a from BAO data alone are 
not good, so the results are not shown and BAO data 
alone is not discussed in the rest of this paper. The con- 
straint from SNe la data alone is similar to that from the 
combination of SNe la and WMAP7 data in the curved 
case. This is easily understood because the addition of 
WMAP7 data in the curved case is used to limit the fife 
around zero. When the BAO data is added to the SNe 
la data, the uncertainty in w a is reduced more than half 
and the FOM becomes 5 times larger. When WMAP7 
data is added to the SNe la data, the uncertainty in w a 
is reduced a little further and the FOM becomes almost 
10 times larger. Compared with the constraint from the 
combination of BAO and WMAP7 data, we see that the 
SNe la data constrains better on wo and the FOM is just 
one half smaller, both BAO and WMAP7 data help re- 
duce the uncertainties in w a , the help from WMAP7 data 
is a little better. When the combined SNe la, BAO and 
WMAP7 data are used, we get better constraints on wq 
and w a . The addition of GRB and H(z) further reduces 
the uncertainties in wq and w a . By using the constraints 
from the combination of all observational data, we re- 
construct w(z) and Om(z) and the results are shown in 
figured ACDM model is consistent with almost all the 
combinations of different data at la level. 

From the above discussion, we find that WMAP7 data 
helps reduce the uncertainties in f2fc and BAO data 
helps reduce the uncertainties in Q m . Neither BAO nor 
WMAP7 data alone gives good constraint on wo and 
w a , but both WMAP7 data and BAO data help SNe 
la data break the degeneracies among the model param- 
eters, hence tighten the constraint on the variation of 
equation of state parameter w a , and WMAP7 data does 
the job a little better. SNLS3 SNe la data alone does not 
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provide good constraints on Q m and Qk, but it provides 
good constrains on the parameters Wq and w a , so it is 
necessary to apply SNe la data to probe the dynamical 
property of dark energy. The addition of H(z) data help 
better constrain the geometry of the universe Qk & n d the 
property of dark energy. Due to the degeneracies among 
the model parameters, we need to measure £l m and f2& 
more precisely in order to better probe the property of 
dark energy. In other words, we need to combine differ- 
ent observational data such as SNe la, BAO, WMAP7 
and H(z) data as long as the tensions among those data 
are not too big. The SNLS3 SNe la data fitting pa- 
rameters a and /3 are consistent at la level for different 
combinations of SNe la with other data. 

3.2. Piecewise parametrization of w{z) 

To see the property of dark energy constrained by cur- 
rent observational data, we apply all the observational 
data outlined in the previous section to the piecewise 
parametrization of w(z) for flat case, 



Vde(z) = (i-n m )(i + z)^ 1+w ^ Y[(i+ Zl 



\3(wi-i—ii)i) 



(15) 

where Zj_i < z < Zi, zq — 0, z\ = 0.1, z 2 = 0.4, 2:3 = 0.7 
and Z4 = 1.4. We also assume that w(z > 1.4) = 
— 1 . Following iHuterer fc Cooravl (|2005f ) , we transform 
the parameters Wi to the de-correlated parameters Wi. 
The results of Wi are shown in the lower right panel 
of Figure |U The res ults are similar with th ose using 
Union2 SNe la d ata (lAmanullah et al.l I2010T) and pre- 
vious BAO data (iGaztanaea. Miguel fc Sanchez! [2009a; 
iPercival et al.l2010D in lGong. Zhu fc Zhul(|2011[ ). and flat 
ACDM model is consistent with the observational data. 

3.3. Modified holographic dark energy model 

In addition to the model independent study, we con- 
sider a specific dark energy model in this subsection. Ap- 
plying the relationship between the mass and the horizon 
of a Schwarzschild black hole in higher dimensions and 
holographic principle, a modified holographic dark en- 
ergy model (MHDE) with Hubbl e hor i zon a s the ultravi- 
olet cutoffjvajyDrop_c^ BothDGP 
model (|Dvali. Gabadadze fc Porratil I2OOOD and ACDM 
model are special cases of this mode l. In this model , 
Friedmann equation is modified as (|Gong fc Lil 120101 : 
IDvali fc Turnenl2003h 

e 2 {z) - (1 - n m - n k - n r )E 5 - N ( z ) 
= o fe (i + zf + n m (i + z f + n r (i + z )\ (16) 

where N is the spatial dimension. So we recover the DGP 
model if N = 4 and ACDM model if iV = 5. This model 
has one more parameter than ACDM model, i.e., there 
are three parameters p = (f2 m , Vtk,N) in this model. 
For the curved case, Qf. ^ 0, again we do not apply SNe 
la, BAO and WMAP7 data alone because the constraint 
is not expected to be good, but consider their combina- 
tions with WMAP7 data. The contours of f2 m -r2fc, flfc-iV, 
and ft m -N are shown in Figure[S][a)-(c), and the la con- 
straints are summarized in table |31 For the constraints 
on Q m and fifc, we see that the combination of BAO and 
WMAP7 data does better than the combination of SNe 



la and WMAP7 data. However, for the constraint on N, 
both combinations get similar results. When we combine 
SNe la, BAO and WMAP7 data, the constraints on the 
model parameters f2 m , VL^ and N are further improved. 
With the addition of GRB and H (z) data, the best fit 
value of fifc is moved toward zero and the upper limit of 
N are reduced a little further. The results also show that 
observational data favors ACDM model more than DGP 
model since larger value of N is favored. The SNLS3 
SNe la data fitting parameters a and (3 are consistent 
for different data combinations. 

For flat MHDE model, we consider the constraints from 
SNe la data, the combination of SNe la and BAO data, 
the combination of SNe la and WMAP7 data, the com- 
bination of BAO and WMAP7 data, the combination of 
SNe la, BAO and WMAP7 data, and the combinations 
of all the observational data, the contours of Vt m -N are 
shown in Figure [5][d) , and the la constraints are sum- 
marized in table |3] By setting fi^ = 0, as expected, the 
constraint on tt m and N from SNe la data alone is simi- 
lar to that from the combination of SNe la and WMAP7 
data in curved case. With the addition of BAO data, the 
constraint on Q m is improved greatly, hencefore improves 
the constraint on N. The effect of the combination of 
WMAP7 and SNe la data is similar to that of the com- 
bination of BAO and SNe la data except that the best fit 
values of the parameters f2 m and N become smaller. The 
effect of the combination of BAO and WMAP7 data is 
similar to that of the combination of SNe la and WMAP7 
data except that the best fit value of f2 m becomes big- 
ger and the best fit value of N becomes smaller. When 
we combine SNe la, BAO and WMAP7 data or all the 
observational data, the results are similar. For all the 
combinations, N > 5 at la level. The SNLS3 SNe la 
data fitting parameters a and j3 are also consistent for 
different data combinations. 

3.4. qi — q2 parametrization 

In this subsection, we reconstruct the deceleration 
parameter q(z) with a simple two-parameter function 
(|Gong fc Wajlal200l . 



*(*) = J 



1 , qiz + q2 



{l + zf 



(17) 



This parametrization recovers the matter dominated 
epoch at high redshift with q(z) = 1/2. The dimen- 
sionless Hubble parameter is 



E{Z)- 



exp [Jq[1 + q(u)]dbx(l + u) 
(1 + z) 3 / 2 exp 



92 1 qiz -q 2 
2 ^ 2(l+z) a 



(18) 



Since E 2 (z) w (1 + z) 3 exp(gi + q 2 ) when z > 1, so 
the role of matter energy density is played by the sum 
of the two parameters, qi + q 2 = In fl m . Although Q m 
and fifc are not model parameters in this parametriza- 
tion, the comoving distance depends on 17^ through the 
function S^, in order to better constrain the model pa- 
rameters p = (91,92), we consider the flat case fife = 
only. As discussed above for the CPL model, the flat 
assumption of Clk — may impose biased prior in the es- 
timation of cosmological par ameters due to the degenera- 
cies a mong Q m , and w (Clarkso n. Cortes fc Bassettl 
l2007f h For this model, the only effect of flk is through Sk, 



() 



and Sk(x) ~ x when fifc is small, so the impact of the flat 
assumption is expected to be small. Fitting this model to 
SNe la data alone, we get the marginalized la constraints 



0.1 



= -1.68 



+0.98 
-0.87' 



<T2 



= -1.06 



+0.19 



a = 1.4271°;$; and 



(3 = 3. 262^0 ioi- The contour plot is shown in Figure 
[6{a). Using these results, we reconstruct q(z) and Om(z) 
and the results are shown in Figures EJb) and [He). So 
q(z) < when z < 0.5 at 2a level and flat ACDM model 
is consistent with the model at la level. With the SNe 
la data alone, the evidence for current acceleration and 
past deceleration is very strong. 

When we fit the model to BAO or WMAP7 data, we 
need to include the radiation-dominated era, and the nui- 
sance parameters fli,h 2 and f2 m /i 2 which are just data fit- 
ting parameters, so we do not apply BAO and WMAP7 
data alone. For approximation, we take the following 
Hubble parameter, 



E 2 (z) = f} r (l+z) 4 +(l+z) 3 exp 



32 



92 



1 



, (19) 



where the current radiation co mponent Q r — 4.1736 x 
10- 5 h- 2 (jKomatsu et all 120111 ). Fitting the model to 
the combined SNe la, GRB, BAO, WMAP7 and H(z) 
data, we get the marginalized la constraints, q\ = 



0.18] 



012 and q 2 

o 9fjo+0.115 
°- zoz -0.103- 



-1.45 ± 0.1, a = 1.423 



+0.119 

-0.096 



and 



/3 = 3.262T£^. The contour plot is shown in Fig. 
[HJa) . Compared this result with that obtained from SNe 
la data alone, we find that they are inconsistent at la 
model, this shows the tension between SNe la, BAO and 
WMAP7 data in fitting this model. Using the qi-q 2 con- 
tour, we reconstruct q(z) and Om(z) and the results are 
shown in Figures EJb) andHJc). We find that q(z) in- 
creases with the redshift and q(z) < when z < 0.5 at 
2a level, flat ACDM model is inconsistent with the model 
at 2a level. These results may suggest that the approxi- 
mation (IT9l) is not good at high redshift. Since fi m does 
not appear in this model, it may not be straightforward 
to apply BAO and WMAP7 data, this needs to be further 
studied. 

3.5. Piecewise parametrization of q(z) 

We also apply the piecewise parametrization to study 
the property of the deceleration parameter q(z). For 
Zi-i < z < Zi, we have 



JY 



E(z) = (l + z) 



1+9N 



IK 1 



Zi-i) 



(20) 



In this model, we have four parameters p = 
(gi, q 2 , 93, 94). We think this model approximates the 
behavior of E(z) in the redshift range z < 1.5. In the 
radiation dominated era, w e add the radiation contribu- 
tion also. Again we follow iHuterer fc Cooravl pOOl to 
transform the correlated parameters qi to uncorrelated 
ones. Fitting the model to all observational data, we re- 
construct the evolution of q{z) and the results are shown 
in Fig. Efd). Similar to that ob tained by Union2 SNe 
la data (|Gong. Zhu fc Zhdl20Tl . we find that q(z) < 
when z < 0.6 and q(z) > at high redshift, so the evi- 
dences for current acceleration and past deceleration are 
very strong. 





Figure 1. The marginalized la and 2a contour plots of Q m and 
Hj, for the ACDM model. The green lines label the constraints 
from SNe la data only, the yellow lines label the constraints from 
BAO data only, the cyan lines label the constraints from the com- 
bination of SNe la and BAO data, the magenta lines label the 
constraints from the combination of SNe la and WMAP7 data, the 
blue lines label the constraints from the combination of WMAP7 
and BAO data, the red lines label the constraints from the combi- 
nation of SNe la, BAO and WMAP7 data, the black lines label the 
constraints from the combination of SNe la, BAO, WMAP7 and 
H(z) data, and the shaded regions label the constraints from the 
combination of all the observational data described in section 2. 




Figure 2. The marginalized la and 2a contour plots for the 
curved CPL model. The meaning of different colors is the same 
as that in Figure [T] The dashed line in the wo-w a contour de- 
notes the condition wq + w a = 0. The + sign denotes the point 
corresponding to the ACDM model. 



4. CONCLUSIONS 

It is well known that BAO data is more sensitive to 
ft m and WMAP7 data is more sensitive to fifc. As more 
data points become available and the data becomes more 
accurate, we are able to constrain the cosmological pa- 
rameters better. The constraints on ACDM model from 
SNLS3 SNe la data alone are similar to that from BAO 
data alone as shown in figure [TJ so we may wonder if 
we can use BAO data alone to constrain the property 
of dark energy. Applying BAO data to CPL model and 
MHDE model, we find that the constraints on the dy- 
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Table 1 

The marginalized lcr errors for U m and fij, in ACDM model constrained by different observational data 



Data 




_2/» 


a 


j3 


SNc la 


17+"' 1 

u - i ' -0.09 


0.15 ±0.25 


1 -* ZJ -(I,094 


3 - 267 -0.109 


BAO 


26 +u uy 

u - ZD -0.03 


-0 16+ U ' ;1S 
u - ±D -o.n 






SNc+BAO 


0.27 ±0.02 


-0 11+ u - uy 
u - 1J --0.07 


1.428+1-- 


^.Z4y_ 095 


SNe+WMAP7 


22 +U - Ub 
u - zz -0.03 


0.01 ±0.01 


i-422l ; ^ 


o 9fi +U.115 
°" 2D -0.106 


BAO+WMAP7 


29+ ,u ' 2 
u - za -o.oi 


-0 004 +u - uuv 






SNe+BAO+WMAP7 


0.28 ±0.01 


-0 004+ lUUI(i 


^•^ -0.088 


9K1+U.1U4 
O.ZD1_q llg 


SNe+BAO+WMAP7+.ff(z) 


n 9S +U.IW 

u - 28 -o.oi 


0.0006±™ 


i 429 +u.iuy 

i.4zy_ () 107 


°-^* J _0 .099 


All 


0.28 ±0.01 


o.ooosi™ 


1 421 +u ' nv 


-} 9-iq + U-J^ 
°- z *°-0.098 



Table 2 

The marginalized lcr constraints on CPL model by different observational data 



Data 




f2 fc 


w 


»« 


a 


P 


FoM 


BAO+WMAP7 


0.35 ±0.07 


-0.023 ±0.013 


0.4 ±1.4 


-8.6+^ 






0.36 


SNe+WMAP7 


n ,0+0.(19 
u-^-o.oe 


-0 03 +(MW 




-3 4+^ i 


1 430 +unl 


3 262+ 0117 
°- zoz -0 .104 


1.67 


SNe+BAO±WMAP7 


u -^ 8 -o.oi 


-0 015 +U - UUv 


-0.8 ±0.2 


-2.02+^6 


1.41U_ 0g2 


„ 245 +U.l3V 
o.zto_ 0g4 


7.29 


All 




-0 004 +l, ' ()()(i 
U.UU4_„ ()()7 


-0.96_£S 


-0.66+?"™ 


1 4349+ uiU9 


3 256 +Uli9 
3,ZOD -0 .102 


9.73 


SNe 


0.31.^ 




-0 8 +U ' 4 
u.o_ () 2 




1 432+ ulua 


o 9 cc+U.1^3 

o.zoo_ 099 


1.15 


SNc+BAO 


n 98+U-U3 
U -^ 8 -0.02 








1 41 o + U.12b 

i - 4i °-0.090 


, 97 ,+u.nu 

o.zio_ 112 


5.94 


SNe+WMAP7 


24+ u - us 

u "^*-0.02 




-0.9 ±0.2 


-i.i±r 4 


1 41 o+U.133 

J-.-iJ.o_Q ggg 


3 - 284 io.l22 


10.39 


BAO+WMAP7 


0.30 ±0.04 




-0.8lt^ 


-0.90ti y 9 5 






2.38 


SNe+BAO±WMAP7 






i -,9+O.^V 
1 - lz -0.07 


0.32_«;& 


1 429+ ulli! 
± -* za -0 .103 


3 261 +U ' iii 

j.zo±_ 109 


12.52 


ALL 


0.28.^ 




-1 02 +() ' IH 


-020 +, f f 5 


1 416 +U - La '' 


3 94fi+ U ' lx! > 
u.Z40_q 09 g 


15.16 



Table 3 

The marginalized lcr constraints on MHDE model by different observational data 



Data 




n k 


N 






BAO+WMAP7 


0.28 ±0.02 


-0.01 ± 0.01 








SNc+WMAP7 


29+ (M,Y 

u - za -0.08 


-0.02 ± 0.03 


6.5+}. 


1.434 ±0.108 


3.268 + 0.11 


SNc+BAO±WMAP7 


u - J8 -o.oi 


-0.011 ± 0.007 




1 415+ u - iM 

J-.-iJ.J_Q ()g8 


O.Z-iO_ () Qgg 


All 


0.28 ±0.01 


-0 001+ U UU4 
u - uul -0.006 


t- n+0.8 
0.3 


1 4i5+U-^« 

— 0.090 


o 9^+0.122 
°- ZJ1 -0.098 


SNe 


0.26 ±0.09 




D - 1 -1.8 


1.435 ± 0.108 


3.268 + 0.11 


SNc+BAO 


U ' 2 °-0.03 




5 8 +i ^ 


, 4 o 2 +U.112 

l.ihz_ 0A04 


„ 254 +U.lii« 
J - !o *-0 .095 


SNc+WMAP7 


25+ (MB 
u - ZJ -0.02 




t- o+u.a 

°- o -0.3 


1 434+ 011 


3.259_of 


BAO+WMAP7 


0.29 + 0.02 




^■ a -0.3 






SNc+BAO+WMAP7 


0.28 ±0.01 




c o+().(i 
o - J -0.3 


1 423+ u - ii8 

1.1ZO_Q Q9g 


3 264 +(, ' il,(; 

o.ZD4_ 114 


All 


0.28 ±0.01 




t- ,+U.S 
o.o_q 2 


1.4±0_Q Q9 


3 9^4+ u ' llb 



namical behavior of dark energy from BAO data alone 
are much worse than those from SNe la data alone. Al- 
though SNe la data alone is not able to provide good 
constraints on f2 m and f^, it provides much better con- 
straint on the equation of state parameter of dark energy 
compared with that from BAO and WMAP7 data alone. 
Since the way that the model parameters are degenerated 
is different for SNe la, BAO and WMAP7 data alone, it 
is a good idea to combine different dataset to get better 
constraint on the property of dark energy. 

For the curved CPL model, we find that the combina- 
tion of BAO and WMAP7 data gives better constraints 
on Vt m and ilk than those from the combination of SNe la 
and WMAP7, but the constraints on wo and w a from the 
combination of BAO and WMAP7 data are much worse 
than those from the combination of SNe la and WMAP7 



data. The variation of the equation of state parameter 
w a is reduced more than half when we add BAO data to 
the combination of SNe la and WMAP7 data, and the 
FOM increases more than 4 times. The variation of the 
equation of state parameter w a is reduced more than 5 
times when we add SNe la data to the combination of 
BAO and WMAP7 data, and the FOM increases more 
than 20 times. Both WMAP7 data and BAO data help 
SNe la data break the degeneracies among the model pa- 
rameters, hence tighten the constraint on the variation of 
equation of state parameter w a , and WMAP7 data does 
the job a little better. GRB data has little effect on the 
constraints when it is combined with SNe la data. H(z) 
data helps move the best fit value of ilk toward zero and 
make the model more compatible with ACDM model. 
For the flat CPL model, we get similar constraints on 




Figure 3. The marginalized la and 2a contour plots of Q m -wo 
and Q m -Wa for the CPL model. The upper panels are for the 
curved case and the lower panels are for the fiat case. The meaning 
of different colors is the same as that in Figure [T] 
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Figure 5. The marginalized la and 2a contour plots for MHDE 
model. The figures from the upper left to lower right are labeled as 
(a)-(d), respectively, (a)-(c) are for the curved MHDE model and 
(d) is for the flat MHDE model. The meaning of different colors is 
the same as that in Figure [T] 
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Figure 4. The marginalized la and 2a constraints from observa- 
tions. In the upper left panel, we show the wo and w a contours 
for the flat CPL model from SNe la data alone (the green lines), 
the combination of SNe la and BAO data (the cyan lines), the 
combination of SNe la and WMAP7 data (the magenta lines), the 
combination of BAO and WMAP7 data (the blue lines), the com- 
bination of SNe la, BAO and WMAP7 data (the red lines), and 
the combination of all data (the shaded regions) . The dashed line 
in the wo-w a contour denotes the condition wo + w a = 0. The + 
sign denotes the point corresponding to the ACDM model. In the 
upper right panel, we reconstruct the evolution of w(z) by using 
the constraints from the combination of all data for CPL model, 
the shaded regions are for flat CPL model, and the black lines 
are for curved CPL model. In the lower left panel, we reconstruct 
Om(z) by using the constraints from the combination of all data 
for flat CPL model, the solid line is obtained by using the best 
fitting values of Q m , wo and w a . In the lower right panel, we 
show the observational constraints on w(z) by using the piecewise 
parametrization. 

wq and w a for SNe la data alone and the combination of 
BAO and WMAP7 data. Replacing the SNe la data by 
the combination of BAO and WMAP7 data, the FOM is 
just doubled. When the BAO data is added to the SNe 
la data, the uncertainty in w a is reduced more than half 
and the FOM becomes 5 times larger. When WMAP7 
data is added to the SNe la data, the uncertainty in w a is 



Figure 6. The marginalized la and 2a constraints on qi-q2 
parametrization and the piecewise parametrization of q(z), the ma- 
genta lines represent the results obtained from SNe la data alone 
and the black lines represent the results obtained from all the ob- 
servational data. The figures from upper left to lower right are 
labeled as (a)-(d), respectively, (a) shows the contour plots for q\ 
and </2i (b) and (c) show the reconstruction of q(z) and Om(z), 
(d) shows the results for the piecewise parametrization of q(z) con- 
strained by all the observational data. 



reduced a little further and the FOM becomes almost 10 
times larger. Both BAO and WMAP7 data help reduce 
the uncertainties in w a , the help from WMAP7 data is 
a little better. ACDM model is consistent with all the 
observational data. This point is further supported by 
the reconstruction of w(z) and Om(z) as shown in Figure 
[4j For the MHDE model, we find that the constraints on 
N are similar for the combination of SNe la and BAO 
data, SNe la and WMAP7 data, and BAO and WMAP7 
data. We also find that ACDM model is favored against 
DGP model. 

To study the acceleration of the expansion of the uni- 
verse, we reconstruct the deceleration parameter q(z) 
with a simple two-parameter function and the piecewise 
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parametrization which approximates the evolution of the 
universe in the redshift z < 1.5. For the SNe la data 
only, we see strong evidence that q(z) < in the red- 
shift z < 0.5. The la contour from SNe la data only is 
inconsistent with that from the combination of all data, 
it seems that there exists some tensions between SNe 
la data and other data, however the inconsistency may 
come from the approximation we made at high redshift 
or the way we apply the BAO and WMAP7 data. Note 
that the usual model parameters O m /i 2 and flbh 2 do not 
appear in the q(z) models, and BAO and WMAP7 data 
depend on those parameters, so we may not apply the 
BAO and WMAP7 data as usual, this needs to be fur- 
ther studied. 

Although BAO and WMAP7 data provide reasonably 
good constraints on £l m and f^, it is not able to con- 
strain the dynamics of dark energy, we need SNe la data 
to probe the property of dark energy, especially the vari- 
ation of the equation of state parameter of dark energy. 
The addition of BAO data helps reduce the error on Q m 
and the addition of CMB data helps reduce the error 
on £l fe , so both BAO and CMB data help SNe la data 
tighten the constraints on the equation of state of dark 
energy due to degeneracies among Q, m , flk and w(z), but 
neither data alone can be used to probe the dynamical 
property of dark energy. For the SNLS SNe la data, the 
nuisance parameters a and /3 are consistent for all differ- 
ent combinations of data. Their impacts on the fitting 
of cosmological parameters are minimal. 
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